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Abstract 



Spin dependent fragmentation functions for heavy flavor quarks to 
fragment into heavy baryons are calculated in a quark-diquark model. 
The production of intermediate spin 1/2 and 3/2 excited states is 
explicity included. The resulting A& production rate and polarization 
at LEP energies are in agreement with experiment. The A c and H c 
functions are also obtained. The spin independent fi(z) is compared 
to data. The integrated values for production rates agree with the 
data. 
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1 Introduction. 

Fragmentation of quarks and gluons into hadrons is an important subject for 
which QCD, in both the soft and hard regions, should be applicable. Be- 
cause of the entanglement of both of these regions, it is difficult to calculate 
fragmentation functions in general. However, when fragmentation involves 
heavy flavor quarks and hadrons the situation is clearer - factorization of 
the soft and hard parts can obtain. Heavy quark physics lends itself to per- 
turbative QCD. The calculations are considerably simplified by the peculiar 
feature of heavy physics just as in atomic physics, where the presence of the 
heavy nucleus, with a mass much larger than the average momentum transfer 
inside of the system, effectively reduces the number of the degrees of free- 
dom. A heavy quark or a heavy atomic nucleus can be considered effectively 
immobile. Of course, such a description is vitiated as soon as the momentum 
transfers become of the order of the heavy mass, but then an expansion in 
powers of the inverse mass works well |]J . 

Until recent years, heavy quark fragmentation functions have been outside 
the scope of experimental verification. Theoretical predictions for fragmenta- 
tion into light flavor baryons have been based on a range of phenomenological 
models 0, [3] and Monte-Carlo simulations M. Subsequently, spin dependent 
fragmentation functions were put on a firm foundation employing light-cone 
field theory || and the theoretical implications of the heavy quark effective 
theory were incorporated 0. As these developments occurred it was realized 
that the masses of the heavy flavor quarks allow for perturbative calculations 
in the case of the doubly heavy mesons |7|, §. Later we extended these calcu- 
lations to include the heavy quark fragmentation of baryons ||, [K|. (At the 
same time a Russian group |11], [12] made a similar extension to baryons. More 
recently the light cone expansion of the fragmentation functions was used by 
the Amsterdam group [I3| with an algebraic model to generate predictions 



similar to ours.) The results of our phenomenological approach were inter- 
esting for several reasons. Working in the rather general parameter space 
we obtained the result that the direct production rate of the spin excited 
baryons is on the order of or higher than that of the ground state baryons of 
the same flavor flt| . We also found a functional dependence quite distinct 
from the simple Peterson parameterization 0]. 

That excited baryon result may turn out to be a sensible way to explain 
the tendency seen throughout the scarce experimental data to indicate the 




Figure 1: Diagram representing the combined perturbative and soft QCD 
calculation of fragmentation through excited baryons. 

high depolarization of the heavy flavor baryons. If most of the singly heavy 
baryons are produced in the ground state, then fragmentation can be seen as 
the heavy quark picking up a couple of light sea quarks. Since the momentum 
transfers associated with the quark fragmentation are not large enough to 
make the heavy quark flip its spin, the expected depolarization is minimal, i.e. 
the spin of the heavy baryon is expected to be nearly the same as that of the 
heavy quark. However, if we include spin excited baryons into consideration, 
the resulting polarization of the final baryon is less clear. The excited states 
will decay strongly into the ground state almost entirely via pion emission 
(with the exception of one radiative decay for the S c (| )) and mix with the 
directly fragmented sample. 

It is our purpose herein to explore the fragmentation into excited heavy 
baryons which then decay and populate the observable ground state fragmen- 
tation. Since we are interested in the spin dependent fragmentation functions 
for ground state heavy baryons, the intermediate state resonance production 
will have a significant effect. In the next section we present the particulars 
of the model calculation, incorporating QCD, heavy quark approximations 
and the quark-diquark model of the baryons. The spin dependent functions 
are developed and discussed in the following section, including the relevant 
baryon wave functions in the quark-diquark model, along with the compar- 
ison with data, where available. The final section discusses the implications 
of these results in a broader context. An Appendix presents the many scalar 
products that arise in the final sum over states. 

2 Calculation details. 

The indirect fragmentation process is shown in the Fig. |I[ The heavy quark 
emits a hard gluon that splits into the diquark pair. Diquarks can be either 
scalar or vector, with the vector diquark having the larger mass (as suggested 
by the nucleon-A mass difference). If the diquark is scalar, the heavy quark 



forms a ground state spin ~ baryon with it. Since the scalar diquark is 
spinless and the ground state has no orbital angular momentum, the heavy 
baryon has nearly the same helicity as the heavy quark. Vector diquarks 
have spin 1 and the resulting helicity of the baryon is not the same as that of 
the single heavy quark, but naturally depends on the helicity of the diquark. 
That does not produce any depolarization of the heavy quark by itself - 
the heavy quark inside of the baryon does not change its helicity. However, 
depolarization may occur during the "relaxation" of the baryon into the 
ground state. The depolarizing features of the spin excited baryon decay into 
the ground state were discussed by Falk and Peskin |f[3| . Two parameters are 



crucial in determining the amount of depolarization relaxation. Let AM be 
the mass splitting between the | and | excited states of the heavy baryon. 
So AM will be related to the spin-dependent part of the QCD interaction 
(which has a scale determined by I/tuq) that, in turn, is responsible for 
flipping the quark spin. The excited baryon states lifetimes will be nearly 
the same, f. The longer this lifetime is, in comparison to the time required 
for heavy quark spin flipping, the more likely it is for the diquark and quark 
to mix together forming a randomized quantum spin state. The time it 
would take for the heavy quark to flip its helicity is proportional to ^g. 
The relative magnitude of these two parameters determines the nature of the 
decay. Three distinct regions in the parameter space can be identified: 

1) r >> AM. In this case the heavy quark would not have enough time 
to "mix" with the diquark or light degrees of freedom and the ground state 
baryon helicity would be the same as the initial helicity of the heavy quark. 

2) T << AM. The decay proceeds very slowly, so the resulting ground 
state baryon is completely depolarized.. 

3)r « AM. This is the case of partial decoherence, the borderline be- 
tween cases 1 and 2. 

It is important to note that the decay width T is independent of the 
heavy quark mass. Indeed, in a simplified model the decay may be due to 
the diquark and light degrees of freedom only. On the other hand the mass 

A 2 

splitting AM oc QCD is inversely proportional to the mass of the heavy 
quark. This shows that in the limit of an infinitely heavy quark mass there 
is no depolarization. However, for any given finite mass of the heavy quark 
all of the above possibilities may occur. 

There is very little experimental data on the absolute magnitude of the 



decay widths for the heavy flavor excited baryons, but theoretical predictions 



based on the heavy quark effective theory and potential models (see Ref. [15 



and the summary of Falk and Peskin JTJ]]) indicate that the decay width is 
of the same order of magnitude as the mass splitting for E& and smaller than 
the mass splitting for E c . That supports baryon relaxation as the source 
for the spin depolarization. We will return to this issue after presenting an 
outline of the calculation. 

We now proceed with the perturbative calculation of the fragmentation 
functions. For simplicity we will first consider the spin independent fragmen- 
tation function fi(z) (we use the notation of ref. || but without the carat 
in /i) . Once the formalism is established and the calculations are made, it 
will be easier to turn to the case of the spin dependent h± and g\ in the next 
section. The partial width for the inclusive decay process Z° — > H + X 
can be written in general for any hadron H as 

dT(Z° -»■ H(E) + X)=Y f dz dt{Z -► UE/z) + X, fj.) D^dz, fi), (1) 

where H is the hadron of energy E and longitudinal momentum fraction z 
relative to the parton i, while /x is the arbitrary renormalization scale whose 
value will be chosen to avoid large logarithms. The function D^h is the 
fragmentation probability in a paricular channel. The parton % is a heavy 
quark and the baryon H = B could be Aq, Eq, or Eq in our model. In order 
to extract the leading twist fragmentation function f\(z,fi) out of the total 
fragmentation probability we take the high energy limit. More precisly the 
limit is of the proper light cone component Iq + I3. As long as we are only 
looking at leading twist functions, the high energy limit produces the same 
result. The expression, Eqn. [I] will be simplified if we restrict ourselves to 
the fragmentation channel shown in Fig. |l]. Then, while z is kept fixed, the 
limit of large mass of the Z° along with large energy of the heavy quark, q , 
and the baryon, l , yields 

lim dT(Z° -> B(E) + X) = lim / dz dt(Z -» Q(E/z) + X, ft) fAz, u). 

lo — >oo go — ►oo./o 

(2) 
and 

/ dzfifa/i) = — , (3) 
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where Ti is the decay width of Z° into the ground state baryon and appro- 
priate remnants - antiquark, spectator diquark and pion, while Tq is the total 
decay width of the Z° into the heavy quark pair |5|. 

The decay width of the (infinitely massive) Z° into the inclusive heavy 
ground state baryon Bq is viewed as a direct decay into Q + Q wherein 
the quark is off-shell and subsequently fragments. The fragmentation is fi- 
nally into the Bq, an anti-diquark D and a single pion (or photon) from the 
intermediate resonant baryon decay and is represented by the integral: 

Tl= 2^I ^ [dl] [dP>] [dW] ( 27r ) V ( Z -q-l-p'-W)^ 2 (4) 

where q, I, p' and 7r are the 4- momenta of the Q, Bq,D and the pion (or 
photon), respectively, and the amplitude Mi is summed and averaged over 
unobserved spins and colors. We use the notation [dp] = d 3 p/(l67r 3 po) fo r 
the invariant phase space element. To isolate the fragmentation function, 
the production of the fragmenting quark (dT of Eq. [I]) must be factored out. 
The fictitious decay width for the Z° — > Q + Q, with the Q-quark on shell is 

To = ^J[dq][d q ]{2^5\Z-q- q )^£W,\ 2 . (5) 

with q the heavy quark 4-momentum. 

In order to factor the fictitious decay width out of Eq. |] we have to 
transform the phase space variables. That can be achieved by introducing 
new, production independent variables x\ = Po ^ PL and %2 = V L that can 
be loosely thought of as Feynman scaling variables for each subprocess, i.e. 
excited baryon production and decay. We will introduce a further simpli- 
fication - the ratio of the narrow decay width to the mass of the excited 
baryons allows the narrow width approximation to be used. The square of 
the denominator for each excited state propagator that enters the squared 
amplitude can be factored out for each decay channel and is proportional to 
( 2_M 2l )+M 2 r 2 " Hence, in the limit of the small decay width to mass ratio, 
that factor can be approximated by jj^5(p 2 — M 2 ), effectively putting the 
excited baryon back on the mass shell. The cross terms that come from the 
cross products of the propagators of the | and | spin channels of the decay 
disappear if the mass difference is much larger than the decay width. 



The resulting phase space integral can be written as: 

Ti = — l — [\dq]\dq](2n)HUZ-q-q) 

1 2M z 256ir 4 J [ yjL m ' V y H > 

M4 



/( \/l T77 \ f 

ds q 9 s q / d(j)d(pdxidx2 \Ai\ (6) 



Here \Ai\ 5(p 2 — M 2 ) = \M-\_\ . The two angles <fi and (p introduced here 
deserve some special attention and are defined carefuly in the Appendix. 
They are associated with the position of the transverse momentum vector 
in two frames of reference. The first is the frame determined by the three- 
momentum of the heavy quark and a fixed vector perpendicular to it, which 
is arbitrary unless it is the quark's transverse spin vector (which enters in 
h\ only). The angle is the azimuthal angle between this plane and the 
transverse momentum vector (relative to the heavy quark direction) of the 
excited baryon. The second plane (or frame if we add the vector perpen- 
dicular to the first two) is constructed out of the three-momentum of the 
excited baryon and the spin vector perpendicular to that three-momentum 
but having no transverse component relative to the first frame. The second 
angle tp is defined as the azimuthal angle between this latter plane and the 
transverse momentum (relative to the excited baryon direction) of the final 
baryon. 

The spin averaged matrix element \Ai\ has no angular dependences, so 
we can safely integrate over the angles. That will not be true for the g\ and 
especially h\. The choice of x\ and X2 helps to keep the integral symmetric 
looking, but unlike the standard scaling variable z = °~|~ L these variables 
are not experimentally observable. Using z = X1X2 we can finally rewrite the 
phase space integral to be: 



.f dsA s,-Mi.^L } f^ dzlAl ^ (7) 

J z 1 — z J x 2 

After factoring out the production decay width we are left with the some- 
what simpler expression for fa 

fi(z,u) — lim / ds o (8) 

J1K ' W 167T 2 90-oo J Sfh X9 Mn 2 ^ > 



The expression above is general for the four body final state. The model 
dependence is hidden inside of the |Al| with the delta function obtained 
in the narrow width approximation integrated out. The final ground state 
baryon can be produced via one of the two intermediate states. If the states 
are separated by a mass gap wider than the decay width they do no not 
interfere with each other. That is the scenario supported by the theoretical 
predictions and the experimental data we have at the moment. It follows 
that the two channels for the indirect baryon production can be considered 
independently. The amplitudes for both of them can be expressed as: 

A -U(l)[ Kl / 2{ *-^ )Al ' 2 | ^3/2^,(1 + M E K /2 1 (g) 

Ti/2 yT^2 J 

where II is the quark production spinor, K is the decay operator for 1/2' (into 
a pion or a photon) or 3/2 baryon, the subscripted A is the corresponding 
production operator for the excited baryon and P^ft + Ms) is the spin sum 
of the 3/2 baryon. 

K1/2 = gi{£>- 7)75 

k$/2 = -92( P -iy 

11 2 

Pp, = -g^ + Tp^lu+T^il^Pv-luP^ + T^pP^u (10) 

where g\ and gi are coupling constants associated with the decays. Their 
exact value is irrelevant in this narrow width approximation for the resonance 
that decays into a single channel - both denominator and numerator of Eq.||] 
contain the second power of the decay coupling constant. 

Note that S c has to be treated differently. The spin excited states are 
widely separated, so the lowest lying spin | S c (2574) state can no longer 
decay via tc. That only leaves the possibility of photon decay with a branching 
ratio of nearly 100%. For the radiative decay we will consider the simple 
electric dipole transition amplitude, g^U{l) /£(j)U(p), where the photon's 
polarization vector 6(7) enters. This leads to a Ky 2 — #3 / replacing the 
value in Eq. [H]. 

The ground state baryon is composed of a heavy quark and a scalar 
diquark. For direct production of the ground state there is one coupling 
constant for the scalar diquark to couple to the gluon field via a color octet 



vector current - a color charge strength, along with a possible form factor 
F s . 

Jt iS) = 9sFs(k 2 )(p + p%S^\£ Se, (11) 

where p and p' are the scalar diquark 4-momenta and k = p' — p. Then the 
amplitude for direct production becomes 

A sl/2 = -^^F s (k 2 )U B g s [k x - 2m d v x ]P x , (12) 



where 



P^ aA% oi m Q (l + v)+k 



The amplitudes for the production of the ~ and | states involve the gluon 
coupling to the vector diquark. For the vector diquark, the color octet current 
(which couples to the gluon field vector) is more complicated [Kj. There are 
three constants - color charge, anomalous chromomagnetic dipole moment 
K, and chromoelectric quadrupole moment A, along with the corresponding 
form factors, F E , F M , and Fq, 

J£ {V) = 9s{\ A ) Pa {F E {e)ie«{p)-S{p')}{p + p'), 

+ (1 + K)F M (k 2 )[e«(p)p ■ efttf) + e^(p')p' ■ e a (p)} 
+^rF Q (k 2 )ie a p ( P yj(p>) + \g pv e\p) ■ e^{p')\k^{p + p% } , 

(14) 
where A is the color octet index, a, j3, ..., are color anti-triplet indices, the 
e's are polarization 4-vectors for the diquarks. The chromoelectric part of 
the matrix element contributing to the spin | baryon is 

Aei/2 = -^=F E (k 2 ) l5 ^^g s e;[k x - 2m d v x ]P x . (15) 

The chromomagnetic contribution to the spin | baryon is taken to be small 
based on earlier extimates from baryon spectroscopy [0, [TO, |TB| . The quadrupole 



coupling is assumed inconsequential due to the more rapid fall-off with mo- 

3 
2 



mentum transfer from dimensional counting rules. For the spin | baryon the 



corresponding chromoelectric amplitude is 

V>(o) 



A E3/2 = -^=F E {k 2 )g s ? v [k x - 2m d v x ]P\ (16) 



After some simplification we can write the above equations including the 
decays: 

U(l)K 1/2 (p + Mz)Ai/2 = j== —F E (k^ 



F(0(i--M s )(M s j? + {e*p)) 

[2MUl - r) - 2^(kp) + M s /c] (17) 



for 7r decay of the 1/2'; 



Tln\K (a^m\a ^(0)(M s + M a )^ 1/2 f 2«? 

^/(Z)^i/2(p + Af E )>li/2 = 7=^ '—F E (k^ 



\?Qm d Mx(s-ml 

U(l)l, K* 2 ){i> - M S )(M S ^* + (e*p)) 



r2/r ... o("P) 



[2M^(l-r)-2^(M + M E ^] 
for the photon (e(7r 2 )) decay of the 1/2'; 

I7(/)^ /2 ^(^ + M s )^ /2 = ^^F^ 2 ^ 



v/2^ ^ v y M E (s-m|) 
Z7(/)(^-P)P^(i> + M s )e^ 



r 2 /1 ., n (np) 



[2M£{1 - r) - 2^(M + M E /;] (19) 

for the 7r decay of the 3/2 state, where r = rrid/Mz is a measure of the 
departure from the heavy quark limit. Note that the H c (| ) radiative decay 
is somewhat different from i^i/2 above, but easily accommodated in the sum 
over intermediate states. 

The fragmentation functions must be independent of the heavy quark 
production mechanism. For simplicity then, we actually calculate the decay 
of a heavy scalar breaking into the heavy quark antiquark pair with the 
quark further fragmenting into the baryon. The remaining calculations are 
straightforward, but quite complex. We will sketch them in the next section. 
The results for the different heavy quark flavors are presented following that. 
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3 Spin dependent fragmentation functions. 

The leading twist fragmentation functions g\ and h\ carry valuable informa- 
tion about the helicity and transversity transfer in the system. The meaning 
of these functions first introduced by Jaffe and Ji || is straightforward. The 
function g\ is the difference between probabilities for having the final baryon 
with helicity aligned with and opposite to the original helicity of the heavy 
quark; h\ similarly represents the difference in the no flip and flip rates, but 
in the sense of the transverse direction of spin (more carefully, the transver- 
sity 17|). 

In the strict limit of the heavy quark effective theory one expects to find 
f\{z) = gi(z) = hi(z) = P5(l — z), with P being an overall production 
rate for the corresponding baryon. This is not quite so for any given finite 
mass of the heavy quark. The only restrictions come from the probabilistic 
interpretation of the fragmentation function (for example, fi{z) > gi{z)) and 
the analog of the proposed structure function inequality |L8| f\(z) + gi(z) > 



\2hi(z)\. Both of these restrictions are satisfied in our model. 

The theoretical prediction made in ref. ]T3|] indicates that the polariza- 



tion of the Aq is heavily dependent on the production rates for the excited 
baryons. We only consider the lowest spin excited baryons £q(| ) and £q(§) 
in the present paper. Other excitations can be included in two ways: con- 
sidering excited diquarks, either radially or orbitally, or considering quark- 
diquark baryon configurations excited radially and/or orbitally. The former 
case creates a lot of theoretical uncertainties, since the diquark is less bound 
and can hardly be considered as a parton. In the latter case, the wave func- 
tions at the origin are expected to be smaller, so that the production rate of 
such states is smaller. In general, such states may provide for corrections to 
our calculations, but the main features of the fragmentation should remain 
unaltered. 

In our model the spin dependent fragmentation functions can be obtained 
directly, using the modified Eq.j|: 



gi(Q,z) = -J—r lim [" 'ds—d^— t^- (20) 



1 „„ r,dx 2 ^ \A Xv +\ 2 -\A ly . 



hi(Q,z) = j lim / ds—d(j)d<p^^ ' y ' (21) 
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with new indices specifying the spin alignment (\y+ >~ |+ > +i\— >)• An- 
gular integration is especially complicated for hi, because matrix elements 
involve spin projections that make them no longer azimuthally symmetric. 
Naturally, the transverse spin vector would set a preferred azimuthal direc- 
tion that manifests itself in the scalar products of the transverse spin vectors 
with other vectors that have transverse components. These scalar products 
are listed in the Appendix. 

The amplitude A is defined in Eq. ^| and in general can be represented as: 

A = UOT1 (22) 

where O is the expression in curly brackets in Eq. |9|. The square of the matrix 
element can then be written as: 

\A la \ 2 = Tr{m la 0^ lG {£> + M A ) 1+7 2 5/ga O) (23) 

with a being the spin projection index and S a being the spin four-vector 
corresponding to that projection. The two cases for this problem are the 
longitudinal or helicity and transverse spin vectors. The longitudinal spin 
vectors are defined unambiguously by defining the four-vector in the rest 
frame of the particle with time component equal to zero and space compo- 
nent equal to the unit vector pointed in the direction of the Lorentz boost 
that would take the particle back into the lab frame. The transverse spin vec- 
tor (transversity ]17|) has to be perpendicular to the direction of the Lorentz 
boost and has zero energy component. It will not be affected by the Lorentz 
boost, so it will still be perpendicular to the three- momentum of the particle 
even in the lab frame. The tip of the spin vector could be anywhere on the 
unit circle that has its origin at the base of the spin vector and is perpen- 
dicular to the three-momentum of the particle. The ambiguity is resolved 
by choosing the positive direction of the transverse spin of the final baryon 
to "align" (or anti-align) with the arbitrarily chosen positive direction of the 
initial quark's transversity. By taking the three-momentum of the quark to 
be along the z-axis and transverse spin direction to be along the x-axis the 
transverse momentum will be defined to be in the xz plane with it's x compo- 
nent positive. Such a definition of the transverse spin vector also naturally 
introduces the new coordinate system with z along the three -momentum 
of the baryon and the x-axis parallel to the transverse spin. In exactly the 
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same fashion we can introduce the transverse spin vector of the intermediate 
baryon. We set up yet another coordinate frame. The angles <fi and if are 
azimuthal angles of the intermediate baryon in the frame of the heavy quark 
and the final baryon in the frame of intermediate baryon (see the Appendix). 

As was mentioned previously, we can proceed in explicitly calculating all 
matrix elements without losing any generality by assuming that the produc- 
tion mechanism is a simple decay of the scalar particle. This way II = v_ a , 
where —a indicates that the spin of the antiquark is opposite of the expected 
spin of the quark. After simplifying the squares of matrix elements and leav- 
ing only the leading terms in the high momentum limit, we end up with 
the scalar products of all the involved four-momenta and spin vectors that 
are given in the Appendix. We next look at the angular composition of the 
resulting integrands. 

In general, the square matrix elements involve scalar products of the 
available four-vectors: 

q,q,p,l,n, s 9 , s«, s l . 

The spin vectors can be either transverse or longitudinal and n is the four- 
vector orthogonal to the quark momentum that enters in the axial gauge. The 
phase space integration variables are 

Xl,X 2 ,(l),<P,Sq. 

We will use the notation p qt throughout this paper, where the first index 
stands for the three-momentum that originated the frame and the second 
index, if there, corresponds to the component of the vector. So p qt is the 
transverse component of the excited baryon three-momentum viewed in the 
first frame. Note that 

p 2 qt = s q xx{\ - x x ) - M|(l - x x ) - m 2 d x x , (24) 

l\ t = Mlx 2 (l - x 2 ) - M\{1 - x 2 ) - m\x 2 . (25) 

In the Appendix we consider each scalar product that can arise in the 
integrands of Eq. EUI and I2T. 



For the f\ there are no spin vectors since we spin average everything. This 
way the only scalar product that will introduce angular dependence is (ql) (as 
shown in the Appendix), the product of the four- momenta of the initial quark 
and baryon. The angular dependence is proportional to cos((p + (f>). This 
term disappears after the single integral over one of the angles - introducing 
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intermediate baryon degrees of freedom breaks the azimuthal symmetry of 
the system, but the integration restores the symmetry. 

The function gi has the same type of angular dependence. The only scalar 
products that produce azimuthal angles are (ql), (s q l) and (siq). All of them 
are proportional to cos(ip + <fi) and the angular dependence is removed by 
integrating over one of the angles, which is not surprising given that the 
azimuthal symmetry should be preserved again. 

In the case of hi the situation becomes more complex. The scalar products 
of the type (s q p) and (sip) are proportional to cos(0) and cos(y) and we still 
have the scalar product (ql) that is proportional to cos(v? + 0). This type of 
angular dependence is removed only after both angles are integrated. The 
system is never azimuthally symmetric - in choosing the arbitrary transverse 
spin direction we break the rotational invariance. 

In order to obtain the final fragmentation function we only have to inte- 
grate over the angles and change all the variables into X\ and X2, which can 
be done in the high momentum limit. This integration is performed numer- 
ically. There are several parameters that need to be specified. The masses 
of the diquarks are taken as (uu, ud) 0.6 GeV/c 2 and (us, ds) 0.9 GeV/c 2 . 
The quark masses are 4.9 GeV/c 2 for the b-quark and 1.6 GeV/c 2 for the 
c-quark. The A& mass is 5.6 GeV/c 2 ; the £& mass is 5.8 GeV/c 2 (the mass 
of the Sfe* can be taken the same, because difference is in the hundredths); 
the A c is 2.3 GeV/c 2 ; the mass of S c is 2.5 GeV/c 2 . The wave functions for 
the formation of the baryons are obtained from the power law potential of 



Eichten and Quigg [19 



V(r) = -8.064 GeV + 6.898 GeV (r x 1 GeV) ai . (26) 

The resulting square moduli of the wavefunctions at the origin, |-^(0)| 2 , are 
0.46 GeV 3 , 0.35 GeV 3 and 0.51 GeV 3 for A fe ,A c and S c , respectively. 

The integrations produce spin-dependent fragmentation functions that 
are defined at the scale /Iq = ttiq + mdi qU ark- To evolve them to higher 
values of the defining scale (or the typical Q 2 ) we utilize the appropriate 
spin-dependent Altarelli-Parisi integro-differential equations as determined 
by Artru and Mekhfi |2(J. 

In Fig. the unevolved fragmentation functions fi,gi,hx are plotted for 
A c as functions of z. The function fi(z) evolved to 5.5 GeV (half of the 
CESR energy) is shown in Fig. [J| along with the CLEO data [Bl| on A c 
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Figure 2: Fragmentation functions for A c at jjlq. At the peak f\ is largest, 
followed by hi and g\. 

production from one of the decay channels, pK~ir + . The data are given for 
bins of x + , which is related to our z, but somewhat different at finite c-quark 
energy. The normalization of the data (which is arbitrary) was adjusted to 
correspond to the normalization of the predicted curve. There is considerable 
scatter about our curve for fi(z). A smooth curve through the data suggests 
a lower z for the peak position. The data are not sufficiently spaced in x + to 
check whether or not there is evidence for the shoulder in our predicted curve. 
Caution is advised in interpreting the CLEO A c production data as a leading 
twist fragmentation function, given that the c-quark is being produced at half 
the T mass. This is far from a large energy compared to the fragmenting A c , 
i.e. tjMa c /My ~ 0.4 which is quite sizeable. 

The 45 GeV spin dependent functions are shown in Fig. |3] for complete- 
ness. Corresponding results for the b-quark are shown for A& in Fig. || and |[ 
The S c fragmentation functions are presented in Figs. 0J^ and ||], for the 
unevolved, 5.5 GeV/c (for comparing with CESR data) and 45 GeV/c (for 
comparing with LEP data). 

It is clear from these figures that the unevolved functions are all sharply 
peaked at high z and get spread out and smoothed out with evolution. The 
double peak structure noted for the direct production of the ground state 
heavy baryons in our previous work |9|, [H| has been moderated by the con- 
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Figure 3: Spin independent function fi(z) evolved to fi= 5.5 GeV. The data 
are from CLEO (21]]. 
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Figure 4: Fragmentation functions for A c evolved to n= 45 GeV. At the peak 
/i is largest, followed by h\ and g\. 
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Figure 5: Fragmentation functions for A& at fi^. At the peak f\ is largest, 
followed by hi and g\. 




Figure 6: Fragmentation functions for Af, evolved to 45 GeV. At the peak j\ 
is largest, followed by hi and g\. 
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Figure 7: Fragmentation functions for S c unevolved. At the peak f\ is largest, 
followed by h\ and g±. 
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Figure 8: Fragmentation functions for S c evolved to 5 GeV/c. At the peak 
/i is largest, followed by hi and g\. 




Figure 9: Fragmentation functions for S c evolved to 45 GeV/c. At the peak 
/i is largest, followed by hi and g\. 



tributions from the excited states. The bound fi + gi > 2\hi\ is satisfied [JT8 
and nearly saturated, as expected in the heavy quark limit. 

The overall rate for a heavy baryon state to be a fragment of the corre- 
sponding heavy quark is obtained by integrating f\(z) over z. The results, 
which include the excited intermediate states are tabulated here along with 
experimental data. 



Total Fragmentation Probabilities 
Particle Experiment 

A c 5.6 ±2.6% [OPAL [22]] 



A b 7.6 ±A.2%[ALEPH 



Prediction 
3.9% 
0.59% 
6.7% 



These results are consistent with experiment within errors. The total frag- 
mentation probability for S c is not given, although the fractional rates for 
different states are known [Q and were stated in our previous work ||10|| . 

The predicted overall net polarization transfer from quark to baryon is 
the integral of gi(z)/fi(z). For these we obtain 0.288,0.239 and 0.302 for 
A c , S c and A&, respectively. The averaged value of the polarization mea- 
surements for A b produced at LEP |25| are —0.45 ± 0.18 compared to the 
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Standard Model expectation of —0.94 for the b-quark polarization. So the 
value of net polarization transfer for b — > A& is 0.48 ± 0.19, consistent with 
our prediction of 0.30. Note that our prediction is based on integrating over 
all z, whereas the experiment necessarily emphasizes a narrower range of z 
centered on the peak production region and there is some fluctuation of the 
value of gi(z)/fi(z) in that region. 

4 Discussion of the results. 

It is known that heavy quark fragmentation into excited charmonium and 
bottomonium states contributes significantly to the high energy hadronic 
production of J/ip, ip' and T states [p6fl . The excited intermediate states 
do not make up completely for the overall theoretical underestimate of the 



cross sections for heavy quarkonium production (color octet [27J schemes 
or alternatives are needed), but they are quite significant. The situation 
in the production of heavy baryons is less clear at this time. However, the 
observation of the depolarization of the heavy quark provides an opportunity 
to pin down the contribution of excited states and to compare the theoretical 
model predictions based on QCD with the experimental data. Such an effect 



for the A& was first estimated by Falk and Peskin [|T4]]. The heavy c and b 
baryons fall into their second or third category. The decay of the excited Ej,, 
E c , 2;, or S c baryons is primarily due to pion emission (with the exception 
of the radiative decay of S c (2574) as we discussed above). The experimental 
data on the decay widths are either unavailable or come with large margins 
of error. They can be theoretically estimated for the heavy baryons [p^] . The 
results clearly indicate that the decay of the excited heavy baryons should 
exhibit at least partial depolarization. 

Such depolarization is more important in the light of the results obtained 
in the direct fragmentation model |5], [CTJ . Since the number of the excited 
states that can contribute is large, they may very well pollute the final sample 
of the ground state baryons and result in an overall depolarization. Based 
on general, model independent reasoning Falk and Peskin |TJ| estimated the 
amount of depolarization: 

MZIJ 2(2 - Wl ) A 

A 6 (+i) 9 + A(5 + 2 Wl ) { ' 
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with all heavy quarks initially taken to be in +~ state. The parameter A 
is the probability of producing the vector diquark (versus the scalar); w\ is 
the probability that the vector is in the helicity +1 or - 1 state. Note that 
when direct production of A b is much bigger than the direct production of the 
excited states, as naively may be assumed to be the case, the depolarization 
is minimal. On the other hand, the signature of our direct quark-diquark 
fragmentation model is the high ratio of the excited state production, which 
leads to a sizeable depolarization. 

The original investigation of excited state contributions by Falk and Pe- 
skin flrjj was made without any dynamical details. What we have done here 



is to go a step further, with an eye toward future high statistics experiments, 
and to look at the dynamics of the indirect fragmentation. We used a per- 
turbative calculation of excited baryon fragmentation and decay. The idea 
was simply to incorporate pion decay into the direct quark-diquark model. 
The basic premise of the model still involves a heavy fast quark originating 
from any high energy source, such as the Z° decay. Following that the quark 
shakes off the gluon that breaks into a diquark antidiquark pair. If the di- 
quark has velocity similar to that of the heavy quark they fragment into the 
baryon. The diquark could be either scalar, and then the baryon produced 
is in the ground state (i.e. A&), or vector, so that the baryon is in one of the 
spin excited states. 

Now note that in Falk and Peskin the final net polarization of Aq (in- 
cluding | and | excited state contributions), obtained from Eq. ^7|is : 

l + (l + 4 Wl )A/9 

p ~ TTa (28) 

That corresponds to the integral of ^- over all values of z. Obviously that 
integrated value is not enough to pinpoint both of the parameters (A and w). 
However, the net polarization of Aq produced purely via decays of excited 
states (direct fragmentation excluded) is: 

P^^l. (29, 

Using this expression to obtain W\ we find: 

For A b : w x = 0.41. Total polarization P = 0.295 
For A c : wl = 0.39. Total polarization P = 0.285 
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In our previous paper we had w± = 0.46 for A& using the direct fragmen- 
tation function of A& and £&. The results here are very close to this value 
although we used different techniques to obtain it. Also, these results assume 
that w and A are independent of z, which generally is not correct (production 
rates may depend on the energies of the gluon and that can be translated 
into z dependence). 

The parameter A has not changed much from the last paper. 

For A b : A = 6. 

For A c : A = 6.3. 

The various fragmentation functions that have been obtained in our model 
have yet to be tested experimentally. The important features of those func- 
tions that should appear in the data include their overall normalizations, the 
resulting longitudianl polarizations, the ratios of the three spin dependent 
functions and the characteristic double hump dependence on z. The func- 
tional dependence on z distinguishes this model from the qualitative param- 
eterization of the Peterson model. The departure from a single sharp peak 
structure is a measure of the departure from the heavy quark limit. So careful 
measurements of the z dependence of any of the functions, f\ (z) , g± (z) or hi (z) 
will be very revealing. 

Appendix 

In this appendix we list all relevant scalar products of two 4- vectors beginning 
with all combinations not involving spin. In each case we show the limiting 
value as the fragmenting quark's energy qo and 3- momentum |lf | become 
large. 

a) (qq) 



(qq) = Qo% + l^f I 2 = y m g + \7tfysq + l^f I 2 + l^fl 2 -»• 2 \~$ I 2 

b) (pq) 

(pq) =Poqo+Pi\~q?\ ->■ 2pi\!f\ = 1x\ |lf| 2 

c) ilq) 

{lq)^2z\!}\ 2 
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d) (nq) 

(nq) =q - \!?\ = \jm\ + |lf | 2 - \~$ 



2|Y 



e) (qp) 

Note that the intermediate quark (momentum q — k) is taken as on-shell 
and the mass of the hadron (momentum p) M-% is approximately m<n quar k + 
m quark- Hence q = k + ^-p and p' = k — j^-. The quantity s q = q 2 . Hence: 

(qp) = rriqM^ + (kp) 



(qq)=m 2 + 2-f(kp) + k 2 = s q 

My, 



P = 2^( M 



(s - m 2 q ) 



So finally 

(qp) = mqM-z + 
f) (nq) 



(Ip) 



(nq) = q + \q\ ->• 2 |lf| 

p = I + W 
(p - I) 2 = m 2 n 

Ml + Ml- 2(pl) = m 2 n 



(Ip) 



M l + M l - m l 



2 

h)(?0 

Evaluating this expression requires careful definitions of coordinate sys- 
tems. Let z = q and x = Sx(q) define the X-Z plane for the incoming quark 
(with transverse spin vector also called s qt ). Let (6,4>) be the polar and az- 
imuthal angles for p. Define a second frame of reference in which z! = p is 
the polar axis and x' = [xcos(9) — z sin(9) cos ((/))} /[l — sin 2 (9) sin 2 (0)] 2 so 
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that x 1 is defined by the excited baryon's transverse spin direction. We can 
express / in the primed coordinates as / = (fo, l P tCOs(ip), l pt sin(ip), l p i). The 
azimuth of q in the primed coordinates is approximately n — <fi. It is more 
precisely = it — + 0(pr ) in the relevant large momentum limit. Using that 
relation we get: 

(~<t • O -»■ \~(?\ T=T1 1 PI + 7=nPqtlpt COS(0' - if) 

\v\ \v\ 

Then after several steps the appropriate limiting value is given by 

{qi) ^ _^ M | + {qp)X2 + M_pAl cos(0 + v) 

Xi Xi X\ 

i) (np) 

(np) =Po+Pqi -* 2xi \lt\ 

J) (nl) 

(nl) = l + l q i -> 1z |lf | 

These complete all scalar products not involving the spin vectors. Now 
we consider transverse spin vectors also. We only have two distinct ones: s qt 
and s lt . 

k) (s qt p) 

(s qt p) = -p qt cosO) 

1) (s qt l) 

This also requires going into the primed frame of p. In that frame: 

(A) = (o, (i + o(^)), o(^), (&£~<« + (i))) 



/j5« , ,s „ / J 



(s 9 */) -> -/ pt cos(y?) - p^jfy cos(0) + O 
(s 9 */) -»■ -/ pt cos(y?) 



p \ \ p 

p qt cos((p) 



X 2 



m) (s qt n) 



,qt r \ . P?*COS(0) 



(s 9 *n) 
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n) (s lt q) 

In the primed frame: 

(A>) = (o, (i + o(-L)), o ( ' ), (-a^ffi + o ( i )/ 

V I p I Ip I p IpI / 

IpI 
1 



o) (s lt q) 



p) (s'*p) 



(A)"" PqtC08(<t>) 

XI 



(s lt q) -»• - p gi cos(0) 

X 1 



(A) ->• Pgi COS(0) 
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